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Abstract. 

We study the entanglement spectrum (ES) of a finite XXZ spin-i chain in the limit 
A — > — 1 + for both open and periodic boundary conditions. At A = —I (ferromagnetic 
point) the model is equivalent to the Heisenberg ferromagnet and its degenerate ground 
state manifold is the SU(2) multiplet with maximal total spin. Any state in this 
so-called "symmetric sector" is an equal weight superposition of all possible spin 
configurations. In the gapless phase at A > —I this property is progressively lost 
as one moves away from the A = — I point. We investigate how the ES obtained 
from the states in this manifold reflects this change, using exact diagonalization and 
Bethe ansatz calculations. We find that in the limit A — > — f + most of the ES levels 
show divergent behavior. Moreover, while at A = — I the ES contains no information 
about the boundaries, for A > — I it depends dramatically on the choice of boundary 
conditions. For both open and periodic boundary conditions the ES exhibits an elegant 
multiplicity structure for which we conjecture a combinatorial formula. We also study 
the entanglement eigenfunctions, i.e. the eigenfunctions of the reduced density matrix. 
We find that the eigenfunctions corresponding to the non diverging levels mimic the 
behavior of the state wavefunction, whereas the others show intriguing polynomial 
structures. Finally we analyze the distribution of the ES levels as the system is detuned 
away from A = — 1. 
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1. Introduction 

Entanglement is one of the most fascinating concepts in quantum mechanics and the 
application of entanglement related ideas to condensed matter in recent years has 
boosted a huge interdisciplinary effort [1]. In particular there has been a growing 
interest in studying the properties of the so called entanglement spectrum (ES) 
[2, 3, 4, 5, 6, 7, 8, 9, 10]. The original motivation for the study of the ES was that 
the ES lies at the heart of the density matrix renormalization group (DMRG) algorithm 
[2]. More recently, the ES has attracted a lot of attention because of its relation to 
low-energy boundary modes in topological phases [3]. 

Considering a bipartition of a system into two parts A and B, the entanglement spectrum 
(ES), is defined in terms of the Schmidt decomposition 

|^) = ^e-^ 2 |^)®|^f). (1) 

i 

where \if>) is the state of the system, and form an orthonormal basis for 

the subsystem A (B). The ES = — logAj} can also be thought of in terms of the 
eigenvalues {A,} of the reduced density matrix pa obtained after tracing out the B part 
of the system density matrix | if> } ( if> | . 

In this work we study the ES of the finite spin-| anisotropic Heisenberg (XXZ) chain in 
the vicinity of the ferromagnetic point A = — 1. (Here A is the anisotropy parameter.) 
In particular we focus on the transformation from the gapless phase at — 1 < A < 1, 
described by a c = 1 conformal field theory (CFT), to the A = — 1 point. In the gapless 
phase (except for the free fermion point A = where many results are available [2]) 
only the distribution of the ES levels is known from conformal invariance [10]. On the 
other hand at A = — 1 the system is not conformal invariant and a dramatic change 
occurs. This is reflected in a peculiar behavior of usual entanglement related quantities 
such as the Renyi entropies and the concurrence [11, 12, 13]. In particular, the Renyi 
entropies in the vicinity of A = —1 show signature of an essential singularity [11, 12]. 
The purpose of this work is to understand how this change in the behavior of the system 
is reflected at the level of the ES. 

The A = — 1 point is unitarily related to the S"?7(2)-invariant Heisenberg ferromagnetic 
chain, and hence the ground state is the degenerate SU(2) multiplet corresponding to 
the highest total spin St- At A = — 1 the states of the multiplet are flat superpositions 
of all the possible spin configurations compatible with total magnetization, S T , i.e., the 
particles are fully delocalized. We therefore refer to the ground state multiplet as the 
symmetric multiplet. In other terms at A = — 1 the ground state shows a "mean field" 
structure, i.e., the system has no notion of distance. This implies for example that there 
is no difference between open and periodic boundary conditions. 

While there are well-known many-body model systems lacking notions of distance, such 
as the Lipkin-Meshkov-Glick model [14, 15, 16, 17] and the Richardson model [18], these 
are generally mean-field like at the level of the Hamiltonian. The specialty of the present 
case is that the XXZ interaction is local, and the symmetric structure appears only at a 
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single value (A = — 1) of the anisotropy parameter. Varying A is a natural way to tune 
away from the symmetric structure. In this Article, we show how this change from the 
symmetric to the more usual situation is manifested in the ES. 

Although the usual interest is in entanglement properties of the ground state, in the 
vicinity of A = —1 it is natural to study the entire set of states which become degenerate 
at that point. The ES properties we report are found to be common to the entire 
symmetric sector. For simplicity of demonstration or explicit calculations, we find it 
convenient to sometimes focus on simpler large- members of this manifold (which in 
the energy spectrum are near the top of this group of states for A > —1) in addition to 
the = ground state. 

At A = — 1, the symmetric structure fully determines the ES [19, 20, 21, 22]. However, 
at A > — 1 the ES is not restricted, and in the limit A — > — 1 + most of the ES levels 
diverge (i.e. most of the eigenvalues of pa vanish). The divergent levels correspond 
to components in the ground state wavefunction at A > —1 that do not reflect the 
symmetric structure present at A = —1. Moreover, while at A = — 1 the ES does 
not contain any information about the boundary conditions or the geometry of the 
partitioning between A and B, for A > — 1 the ES for periodic and open boundary 
conditions display striking differences. For open boundary conditions the ES is less 
dense (some of the ES levels are missing). This is related to there being two boundaries 
between A and B for periodic boundary conditions (and only one for open). This 
difference is also present in the behavior of the entanglement entropy in the gapless 
phase (as 1/3 log I for periodic and as 1/6 log I for open boundary conditions with £ the 
length of subsystem A) . We also find that the multiplicity count of the ES in the vicinity 
of A = — 1 (numbers of ES levels diverging at different rates) seems to be described by 
an elegant combinatorial formula which we present as a conjecture. 
We also analyze the entanglement eigenf unctions, i.e. the eigenstates of pA- While the 
eigenf unctions corresponding to the lower (non-diverging) ES levels mimic the behavior 
of the ground state wavefunction (i.e. they exhibit a symmetric structure in the limit 
A — > — 1 + ), the others show richer structures and in some cases their leading behavior 
(in e = A + 1) can be given in terms of known polynomial functions. 
Finally we investigate the distribution of the ES levels in the vicinity of the ferromagnetic 
point. Although in the thermodynamic limit this is expected to be a universal function, 
for finite chains the ferromagnetic point introduces significant corrections in the whole 
region A < 0. 

For our analysis we take advantage of the fact that the XXZ is an integrable model and 
all its eigenstates can be in principle calculated using Bethe ansatz. We also use exact 
numerical diagonalization. Since ES levels are determined often up to values of ~ 100, 
the reduced density matrix eigenvalues Aj = need to be calculated with very high 
precision. Therefore, as in Ref. [4], we use arbitrary precision numerics. 
In section 2 we discuss the spin-| XXZ chain around A = — 1, showing the ground state 
manifold and how the wavefunctions of this symmetric sector behave at A — > — 1 + . In 
section 3 we present the main features of the ES, reviewing the exact structure known 
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at A = — 1 and then presenting numerical diagonalization data for A = — 1 + e. In 
section 4 we treat the ES perturbatively in the distance A + 1 from the ferromagnetic 
point, using the Bethe ansatz. The combinatorial aspects of the ES are also addressed 
and characterized. Section 5 is devoted to the study of the entanglement eigenfunctions. 
In section 6 we investigate how the distribution of the ES levels changes as the system 
is detuned away from the ferromagnetic point. Details of the Bethe ansatz calculation 
appear in the Appendix. 

2. The XXZ in the limit A — > — 1 + : the symmetric multiplet, energy 
spectrum and wavefunctions 

The spin- 1 XXZ open chain of L sites is defined by the Hamiltonian 

1 L-l L-l 

u = 2E(& + W + A E 5 ^ + i ( 2 ) 

i=l i=l 

and the periodic chain is obtained by adding terms connecting sites L and 1. At A = — 1 
the model is related to the ferromagnetic Heisenberg model (— ^ Si ■ S i+ i) by the local 
unitary transformation S x,v — > —S x,y applied on every second site of the chain |. The 
ground state of the Heisenberg ferromagnet is in the SU(2) multiplet with highest 
total spin St (symmetric multiplet), i.e. the one containing the fully polarized state 
\F) = |ttt ' ' ')• The L + l states of the degenerate multiplet are obtained from \F) by 
successive applications of the lowering operator S? = ^ S~ . It is conventional to refer 
to \F) as the "vacuum" and the overturned spins as 'particles'. The states (S T ~) M \F) 
are equal weight superpositions of all the possible spin configurations with M down 
spins (particles). The magnetization is S T = L/2 — M in terms of the particle content 
M; either (S T or M) is a good quantum number for all A and labels the states of the 
multiplet. 

In Figure 1 (a) we show the energies of the states which are connected to the symmetric 
multiplet at A = —1, as a function of A + 1. The energy levels in the vicinity of A = — 1 

are 



L 1 (L \(L 

S T )(^ + St 



.4 L-1V2 ~VV2 - ( A + 1 >- < 3 > 

The ground state at A > —1 is the one with S T = 0, while the other states of the 
multiplet are excited states. In Figure 1 (b) we plot the wavefunction of the eigenstate 
with two particles (which is part of the symmetric multiplet but not the ground state) 
showing wavefunction components versus the relative position [x<i — X\) of the two 
particles. While at the point A = — 1 the wavefunction is completely flat (the two 
particles are completely delocalized), for A > — 1 this property is progressively lost 
farther away from A = — 1. 

% Since this transformation is a product of local unitary transformations, it docs not affect the 
entanglement spectrum. 
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Figure 1. (a) Energy spectrum around A = — 1 of the XXZ spin chain with L = 6 
and periodic boundary conditions. We show the energies of all the chain eigenstates 
within the symmetric multiplet and > 0. The energy at A = — 1 is E = —L/A. 
The dashed lines are from Eq. (3). (b) The wavefunction of the state with only two 
particles (M = 2, $p = L/2 — M) for a chain with L = 24 and several values of 
A. We plot the wavefunction components versus the interparticle distance X2 — x± 
(xi > x\). The points are exact diagonalization data while the dashed line is the 
Bcthe ansatz result, Eq. (A. 4)). The dashed-dotted line is the asymptotic (A = — 1) 
value, v /2/[L(L - 1)]. (c) The wavefunction of the state with M = L/2 (L = 12) 
particles (ground state at A > —1) plotted versus the relative distance of the first two 
particles. Data are obtained using Eq. (A. 5). We show the wavefunction for three 
different values of A + 1, namely A + 1 = 5 • 10~ 2 , 10~ 2 , 10~ 3 . At each fixed X2 — x\ the 
different points correspond to the possible configurations of the other L/2 — 2 particles. 



This behavior is true for any state in the symmetric multiplet, although more 
complicated to display. For example in Fig. 1 (c) we show wavefunction components of 
the ground state (<Sf. = 0) for several values of A. Now there are multiple components 
for each value of X2 — X\. It is clear, however, that in the limit A — > — 1 the wavefunction 
becomes flatter and eventually becomes an equal weight superposition. 

3. Entanglement spectrum (ES) of the symmetric multiplet: general 
features 

In this section we describe how the structure of the ground state in the vicinity of A = — 1 
determines the behavior of the ES. To calculate the ES we consider a bipartition of the 



Entanglement spectrum of the Heisenberg XXZ chain near the ferromagnetic point 6 



chain into two parts A and B (of lengths L A and Lb respectively). The wavefunction 
\ip) of the chain, given two orthonormal bases for the subsystems \<pf), \ ff) can be 
written as 

hi 

The matrix M^- is a central object in the study of ES. The reduced density matrix 
eigenvalues \ = e~^ [Eq. (1)] are the square of the Schmidt (singular) values of M. 
Since the block magnetization S A is a good quantum number for the matrix M we can 
label the ES levels with S A . 



3.1. The points 



At A = — 1 the flat structure of the states in the symmetric sector, makes it possible 
to construct the matrix M explicitly and to calculate its singular values Any 
eigenstate in the symmetric multiplet is identified given the length of the chain L and 
the magnetization Sf, (particle content). For each of these eigenstates there is one ES 
level per block magnetization sector S A which is given by [19] 



AL,L A } 



} 



log 



L—Lj 



( 



(5) 



As a consequence of the delocalized structure of the states in the symmetric multiplet, 
Eq. (5) holds for both periodic and open boundary conditions. 

After restricting to Sf. = 0, that is the ground state at A > —1, and taking the 
thermodynamic limit L — > oo one obtains 



ALa} 



-log 



L 



La 
2 



(6) 



which has also been obtained in [21]. 



3.2. The ES in the limit A — > — 1 + : exact diagonalization results 

In Fig. 2 (a) we show the ES in the vicinity of A = — 1 (i.e. A + 1 = e = 0) for a block 
with L A = 6 in a chain with periodic boundary conditions and L — 18 sites. For clarity 
we show only the sector of the ES with S Z A = 0, since the result in other S A sectors is 
qualitatively similar. 

The ES levels show singular behavior. In the limit A — > — 1 + there is only one level 
converging to a finite value given by Eq. (5), while all the others diverge. The diverging 
behavior in the ES signals that in the limit A = — 1 + the wavefunction has to be a 
symmetric superposition of all allowed configurations. The components in the total 
wavefunction which reflect this symmetry give rise to a non diverging level, while the 
others contribute to the diverging ones. In Section 5 we will show how this distinction is 




Figure 2. Entanglement spectrum from exact diagonalization (L = 18, La = 6), 
plotted against A. We show only the S A = sector. In the limit A — >• — 1 the lowest 
ES level (in black) converges to the finite value given by (5), while all the other levels 
diverge, (a) Periodic boundary conditions; wide range of A. (b) Periodic boundary 
conditions; zooming around A = — 1. Log-scale on horizontal axis highlights the 
logarithmic divergence as —a log(A + 1) +/5 of the levels. There are multiple ES levels 
for each value of a. We use the notation a[m per ] where m per is the multiplicity of the 
set of levels with same exponent a. (c) The same as in (b) but for a chain with open 
boundary conditions. 



manifested in the structure of entanglement eigenfunctions corresponding to diverging 
and non-diverging ES levels. 

In Fig. 2 (b) we demonstrate that the behavior of the levels in the vicinity of A = —1 is 
given by —a log (A + 1) + /3, a being an integer. The slopes a, going from bottom to top 
in the ES, are given by the sequence 0, 2, 4, 6, ... where refers of course to the lowest 
level. A diverging level in the ES as — alog(A + 1) corresponds to an eigenvalue of the 
reduced density matrix which vanishes as (A + l) a . Focusing on the number of levels 
diverging with the same slope, an intriguing multiplicity sequence appears as higher 
levels in the ES are considered. In Fig. 2 (b) we use the notation a[m per }, m per being 
the number of levels with slope a. We observed that m per (a) is given by the sequence 
m per (a) = 1, 1,2,3, . . .. 

In Fig. 2 (c) we show the same plot but for a chain with open boundary conditions 
(L = 17). The ES shows the same qualitative behavior, although now the slopes are 
a = 0,4,8,12.... The multiplicity structure apparently is described by the sequence 
m open {a) = m per (a/2). 
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Figure 3. ES for the member of the symmetric multiplet containing only two 
particles (M = 2), plotted versus A + l for periodic (PBC) and open (OBC) boundary 
conditions. The block size is La = 6. S Z A (Ma) is the magnetization (particle content) 
of block A. The behavior of the diverging ES levels is £ = — alog(A + 1) + j3. The 
dashed lines are guides to the eye. 



In a generic sector with S Z A ^ 0, the slope sequence (the a's) remains the same 
(a = 0, 2, 4, . . . for periodic boundary conditions and a = 0, 4, 8, . . . for open boundary 
conditions), but the multiplicity structure is different. 

The outlined scenario depends only on the symmetric structure of the states in the 
limit A — > — 1 + , thus the ES obtained from any state within the symmetric multiplet 
is expected to show similar features. This observation is supported in Fig. 3 where 
we show the full ES obtained from the state in the symmetric multiplet with only two 
particles (which is not the ground state for A > — 1, see Fig. 1). We show the ES for 
both periodic (left) and open boundary conditions (right). In both cases there are only 
three possible sectors S Z A determined by the number of ways of distributing the two 
particles in A and B. Nonetheless, the main features observed in Fig. 2 are present: 
(a) for each value of S Z A there is one ES level converging to a finite value in the limit 
A — > — 1 + and this value is given by Eq. (5); (b) the slope of the only divergent level 
is 2 for periodic boundary conditions and 4 for open, in agreement with the sequences 
shown in Fig. 2. 

4. The ES at A — > — 1 + : perturbative and combinatorial structure 

The structure of the ES outlined in the previous sections can be quantitatively 
understood by expanding the state wavefunction for small e = A + l and then calculating 
the ES. For the XXZ we can exploit the fact that the model is integrable and use the 
exact expression of the eigenstate wavefunctions (see Appendix A for details). 
For this description, we choose the orthonormal basis = {\xi, %2, ■ ■ • , xm a )} 
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for the A block, where \xi,x 2 , ■ ■ ■ ,xm a ) are the block configurations with particles 
at positions Xi,x 2 , ■ ■ ■ ,xm a - Similarly for the B subsystem we have {|y?f)} = 

{\ X M A +1, Xm a +2, • • • , Xm }}• 

With this notation the elements of the M matrix [Eq. (4)] are 

M [xm a+ i,-/ m } = x 2, ■ ■ • , x Ma ,x Ma+ i, ■ • • , x M ) (7) 

where [x\, X2, ■ ■ ■ , Xm a ] an d [xm a +i, • • ■ , xm] act as row or column indices. Here 
a(xi, X2, ■ ■ ■ , Xm a , Xm a +i, ■ • ■ , Xm) is the amplitude of the chain configuration with 
particles at positions x±, X2, • • ■ ,Xm a ,Xm a +i, ■ ■ • , Xm- It turns out that the first order 
in A + 1 for the amplitude a in (7) can be obtained analytically for arbitrary M and 
system sizes L and the resulting ES calculated. 

In the three subsections below, we first treat the simplest nontrivial member of the 
symmetric multiplet (subsection 4.1), namely the one with only two particles (M = 2), 
and provide the first order expansion in e = A + 1. In 4.2, we discuss M > 2 and higher 
orders in e, providing a conjecture for the multiplicity sequence of the ES at different 
orders. In 4.3 we comment on the open-boundary case. 

4-1. The sector with two particles (M = 2): first order in A + 1 

In the sector with two particles in the full chain wavefunction, using the expansion 
(A. 4) and (7) the M matrix can be written explicitly. The possible values of Ma are 
Ma = 0, 1, 2. Moreover the blocks of M corresponding to Ma = 0, 2 are lxl (i.e. just 
numbers). In Fig. 3 (PBC) these would give the two ES levels with S A = 1, 3, which are 
not diverging in the limit A — > — 1 + . For these two levels the first order contribution to 
the wavefunction gives a first order renormalization of the result (5). After restricting 
to the case La = Lb = L/2 to simplify the expressions that would be otherwise very 
cumbersome, the ES level in the sector with Ma = 2 is written as 



2 log 



1 jL-2(, L(L-2)_ {A 



(8) 



2 V L - 1 V 24(L - 
Since we chose La = L B the ES level with Ma = is the same. 

More interesting is the sector with Ma = 1 (S A = 2 in Fig. 3), in which case there is a 
divergent level. The block in matrix M describing this sector reads 



M 



L(L-l) 



(9) 

L — 1 L 1 



where ® is the Kronecker product and we defined the three vectors vj^ = k l and 
e = A + 1. The factor A 2 ensures the normalization of the full chain wavefunction 
(see Appendix A). M is a symmetric matrix and quadratic polynomial as a function of 
the coordinates of the two particles. Since the vectors i«w are linearly independent, the 
rank of M is three. However, while two of the singular values of (9) are 0(1) and 0(e), 
the other one is 0(e 2 ) and it is not meaningful because the wavefunction is accurate 
only up to first order. The third singular value is thus "fake" and the effective rank 
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of (9) is two. (See Ref. [4] for a detailed discussion of such spurious eigenvalues in 
perturbative ES calculations.) One ES level of (9) is given by 



2 log 



^(L + 2)(L-2) 

12 



(10) 



24v/2 (L(L-l)y. 

which diverges in the limit A — > — 1 and corresponds to the divergent ES level in Fig. 3 
(PBC). The non diverging ES level in this sector up to first order is given by 



2 log 



(L-2)[L(L-l)Jf (A 



2L(L-1) 12V2(L-1) 

Note that the correction to the singular values, at first order in e = (A + 1), is always 
of type ~Le for large L [Eqs. (8), (10), (11)]. Thus Le = L(A + 1) needs to be small for 
these expansions to be meaningful. 

4-2. Many particles (M > 2) and higher orders in A + 1: combinatorial structure of 
the ES 

In this section we discuss how the structure of the ES for any state within the symmetric 
multiplet builds up as the particle number (M) of the wavefunction increases and higher 
orders in A + 1 are taken into account. 

In each Ma sector the ES levels are organized in powers of A + 1: the lowest ES level 
converges to a finite value in the limit A — > — 1, while the higher levels diverge as 
—a log(A + 1) + (3 with a = 2, 4, 6, .... as observed in Fig. 2. 

Within a given S Z A sector and for each value of the slope a there is more than one level, 
i.e. a non trivial multiplicity structure arises as higher levels in the ES are considered. 
It is natural to investigate how the multiplicity sequence depends on the particle content 
(M). To this purpose in Fig. 4 we provide examples of multiplicity sequences in the 
ES for Ma = 2 and M = 3,4,5 (we consider a chain with L — 18 and La = 8). We 
note that the sequences can be obtained as particular cases from the expansion of the 
so called q-deformed binomial which is defined as: 

r—s+l 

C q (r,s)= ([) = -*^ (12) 



n(l-g fc ) 

k=l 



For example, if we fix M = r = 5, Ma = s = 2 we obtain C q (5, 2) = 1 + q + 2q 2 + 
2q 3 + 2q A + q 5 + g 6 . The sequence 1, 1, 2, 2, 2, 1, 1 given by the coefficients of the different 
monomials matches the one shown in Fig. 4 (c). It can be checked that C 9 (4, 2), C q (3, 2) 
give the other sequences in Fig. 4. We also checked that the result holds for other values 
of M, Ma- Thus we conjecture that for all the states within the symmetric multiplet 
the multiplicity of the ES level diverging as — alog(A + 1) is given by the coefficient of 
q a l 2 in the expansion of the q-deformed binomial (12) with r = M, s = Ma- This result 
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Figure 4. ES for a chain of L = 18 sites, La — 8, periodic boundary conditions and 
several values of A (arbitrary precision exact diagonalization data). M (Ma) is the 
number of particles in the full chain (A block). Logarithmic scale on the horizontal 
axis highlights the divergent behavior (£ = —a log(A + 1) + /3 with a = 0, 2, 4, ...). We 
use the same notation a[m per ] as in Fig. 2 for divergence exponent and multiplicity. 



holds provided that M < La §• We observe that for La < M < Lb the multiplicity 
sequence is given by (12) with r = La, s = Ma, while if M > Lb it is r = La — (M — Lb) 
s = M a -(M-L b ). 



4-3. Open boundary conditions 

The XXZ chain with open boundary conditions can also be treated using the Bethe 
ansatz (details in Appendix A). The first order expansion of the matrix M is 

M [x Ma '+i,-,xm\ = 1 + e[Mi{xx,xz, x Ma ) + M 2 {x Ma +i, • " , ^m)J (13) 

where M.\ and M.2 are polynomial functions that can be obtained from (A. 8). We stress 
that (13) does not contain any cross terms of the form XiXj with Xi G A and Xj £ B 
(note the difference with the result (9) for periodic boundary conditions). This means 
that at first order in e there is no entanglement between A and B. 
Despite the different values of a the multiplicity structure in the ES shows similarities 
with the case of periodic boundary conditions. In Fig. 5 we show the ES levels in the 
sector with Ma = 2 and total number of particles M = 3,4, 5. As already observed in 
Fig. 2 all the ES levels with a = 2k and k odd are missing. Nonetheless the multiplicity 
sequence m open (a) is simply related to the one observed for periodic boundary conditions 
as m open (a) = m per (a/2) with m per given in terms of (12). 

§ Here we arc considering the case La < Lb otherwise the condition is M < min(LA, Lb). However, 
since the reduced density matrices of the two subsystems are isospectral we can always choose subsystem 
A such that La < Lb- 



Entanglement spectrum of the Heisenberg XXZ chain near the ferromagnetic point 12 



M=3 M A =2 



50 



o 

i 25 



[1], 



4[1], 



(a) 



10 



10" 10"' 

A+l 



M=4 M =2 

A 



100 



50 



(b) 



16[l]. s 

N 

12 [1]*, 
8[2] f 

4[1].^._ 



. [l]«-ft »■« ft 

mI i * i * i * 



o 



M=5 M =2 

A 



150 r 



100 



50 



24 [1] . 

20 [1] ft 
16 [2] | 

12 [2] 

4[1] _ ^ 

. oti] . _t 



(c) 



10 



10 



10 10" 

A+l 



10 



10"' 

A+l 



o 
ft 
9 

8 



10" 



Figure 5. ES for a chain of L = 18 sites, La = 8, periodic boundary conditions and 
several values of A (arbitrary precision exact diagonalization data). M (Ma) is the 
number of particles in the full chain (A block). Logarithmic scale on the horizontal 
axis highlights the divergent behavior (£ = — alog(A + 1) + /3 with a = 0, 4, 8, ...). We 
use the same notation a[m per ] as in previous figures 2 and 4 for divergence exponent 
and multiplicity. 



5. The entanglement eigenfunctions at A — > — 1 + 

In this section we examine the entanglement eigenfunctions (eigenvectors of reduced 
density matrix). We show contrasting behavior for the eigenfunctions corresponding to 
the ES levels which diverge in the A — > — 1 + limit and those that do not diverge. 
We discuss the case of a chain with periodic boundary conditions. We provide 
information about structure of the entanglement eigenfunctions in two ways. First, 
up to first order in A + 1 the Bethe ansatz is tractable for calculating the entanglement 
eigenfunctions for arbitrary M, Ma, L (Section 5.1). Second, in the Ma = 1 sector, 
numerical data reveals a polynomial sequence for the forms of the wavef unctions, which 
we express in terms of the discrete Chebyshev polynomials (Section 5.2). 

5.1. First order in A + l 

We start discussing the sector with M — 2. In this sector the entanglement 
eigenfunctions reflect the polynomial structure of the matrix M (see (9)). 
In Fig. 6 (a) (b) we show the ES eigenfunctions for M = 2 and Ma = 1 obtained from 
the expansion (9). The lowest ES level eigenfunction ((a) in Fig. 6) has a parabolic 
dependence on X\ (the position of the particle inside block A) and as the ferromagnetic 
point is approached it becomes flat. The non-diverging ES levels thus mimic the 
structure of the full chain eigenstate. The same behavior is present in the sectors 
with Ma = 0, 2, where the ES has only one (non diverging) level. In contrast, the 
eigenfunction of the first diverging ES level shows linear behavior as a function of X\ 
(Fig. 6 (b)). 
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Figure 6. (a,b) Entanglement eigenfunctions in the sector with only one particle in 
the block A (Ma = 1) and two particles in the full chain wavefunction (M = 2) 
at different values of A. On the x-axis x\ is the position of the particle inside 
block A. (a) is the eigenfunction corresponding to the lowest non-diverging ES level 
while in (b) we show the eigenfunction for the first diverging ES level (see Fig. 3). 
(c,d,e,f) Entanglement eigenfunctions in the sector with two particles in the block 
(Ma = 2). In (c)(e) we show the entanglement eigenfunction for the lowest ES level 
versus Xcm = (#i + x n)l^ an d x r = x 2 — X\ where X\ x 2 arc the position of the two 
particles. In (d)(f) we show the entanglement eigenfunction of the first divergent ES 
level. 

In Fig. 6 (c)-(f) we show the entanglement eigenfunctions for the case with two particles 
in the block A (Ma = 2) and only one in B (M = 3). In (c) and (d) we plot the 
eigenfunction components versus the "center of mass" coordinate (x\ + x 2 )/2 of the two 
particles in A. The entanglement eigenfunction of the lowest ES level (c) exhibits the 
flattening behavior characteristic of the chain wavefunction, while the eigenfunction in 
(d) is again linear. The same qualitative behavior is observed in the relative distance 
X2 — x\ (see Fig. 6 (e) (f)). 

This result is general for all Ma, M sectors — the entanglement eigenfunctions for the 
non-diverging ES levels flatten in the limit A — > — 1 + , while the eigenfunctions for the 
ES diverging as — log (A + 1) are linear functions of the coordinates x±, x 2 , ■ ■ ■ , xm a - 
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Figure 7. ES and entanglement eigenfunctions in the sector with Ma = 1 (exact 
diagonalization data), (left) The ES versus A + 1 for a block with La = 12 embedded 
in a chain of size L = 24. We show the ES levels in the sector with only one particle in 
the block (Ma = 1) and five particles in the full chain (M = 5). The lines are fits to 
— alog(A + 1) + (3, (right) For each level in the ES we show with the same color the 
components in the corresponding entanglement eigenfunction for A + 1 = 10~ 3 . Here 
x is the position of the particle in block A. The lines are the polynomial behaviors 
given by Eq. (16). 



5.2. Entanglement eigenfunctions in the sector with Ma = 1 and arbitrary M 

In the Ma = 1 sector the number of ES levels is given by the total number of particles 
M [Eq. 12]. These levels are shown in Fig. 7 (left) for Ma = 1 and M = 5. For each 
of the ES levels we report (right) the components of the corresponding entanglement 
eigenfunction versus the position of the particle in block A. 

We observe that in the limit A = — 1 + the entanglement eigenfunction of the level 
appearing at order n (i.e. diverging as — alog(A + 1) with a = 2n) is a polynomial 
of order n (as a function of the position of the particle in block A). This is the only 
information needed in order to know the exact functional form of the entanglement 
eigenfunctions. The polynomial can be fixed by imposing that the eigenfunctions are 
mutually orthogonal and are normalized. 

For instance the entanglement eigenfunction of the lowest ES level is the flat 
superposition \ip ) = J2 x c\x), with \x) the configuration with the particle at position 
x in block A. The coefficient c is determined by imposing the normalization of the 
eigenfunction. For the first diverging ES level in Fig. 7 (the red-dotted line) one has 
\ipi) = J2 X ( C ° + c 1 ^)!^). The coefficients c ,^ are fixed by imposing that (ip \ipi) = 
and that \ipi) is normalized. This gives 
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The next diverging ES level (the black-dashed line) is instead 

|V?2) = \l L A (L\-5L\T T) X ^ ( 2 +^+3^A+6a; 2 -6(L A -l)x) |x)(15) 

The procedure allows to obtain the leading order of all the entanglement eigenf unctions 
in the sector Ma = 1. The set of polynomials obtained are known as the discrete 
Chebyshev polynomials T n [23] and are defined as 

r .(«)=^ g) (16) 
where V is such that T>f(x) = f(x + 1) — f(x). 

6. Away from A = — 1: the ES distribution in the gapless phase (A > —1) 




Figure 8. Distribution of the ES for a block with La — 8 embedded in a chain 
of L = 18 (M = 9) sites and several values of A (exact diagonalization). We denote 
with n(£) the number of ES levels smaller than £. The line is the asymptotic CFT 
result (17). 



In this section we show how the symmetric structure of the A = — 1 point affects the 
distribution of ES levels for A > —1. We use periodic boundary conditions and consider 
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M = L/2 (i.e. the state in the multiplet that is the ground state at A > — 1). It has 
been shown [10] that for ID gapless systems described by a conformal field theory (CFT) 
the asymptotic distribution of the ES levels is given by 

n(0 = WVKZ-Un)) (17) 

where n(£) is the number of ES levels smaller than £, £ min is the value of the lowest ES 
level, b = c/61og [L sm(7r / LL a)], c is the central charge (c = 1 for the XXZ spin chain) 
and Iq the modified Bessel function of order zero. This ES distribution depends only on 
the central charge (in this sense is super- universal) , i.e. (17) holds irrespectively of the 
value of A in the gapless phase. 

However, the ferromagnetic point A = — 1 (where the model is not conformal invariant) 
is the endpoint of the c = 1 conformal line which describes the XXZ at — 1 < A < 1. 
Although in the thermodynamic limit (when all the length scales are sent to infinity) the 
CFT result is expected to describe the ES distribution in the whole region — 1 < A < 1, 
for a finite chain deviations from the CFT prediction are expected especially in the 
vicinity of A = — 1. At the A = — 1 point itself, the von Neumann and the Renyi 
entropies show scale invariance but not conformal invariance [21]. The asymptotic 
behaviors are [21] 

W - \ > og m + + om A) , s A = i > og (^) + 1 + om A) 

(18) 

Although the scaling is logarithmic, the prefactors differ from the CFT prediction that 
would give S { 2 ] = f (1 + i) log Lj^ + 0(1) and S A = § log La + 0(1). 
It is interesting to check how this is reflected at the level of the ES distribution. In 
Fig. 8 we show the ES distribution for a chain of L = 18 and block size La = 8. Finite 
size effects are present, but we still notice a good agreement between the numerical 
data and the CFT prediction (17) for A > —0.15. The numerical data in the range 
—0.15 < A < 1 collapse on the same curve, which is well described by the CFT 
result. However, in the region A < —0.15 no data collapse is observed and the CFT 
prediction is a quite poor approximation. Thus for these block sizes and system sizes 
the ferromagnetic point A = — 1 appears to dominate the physics of the whole region 
A < 0. 

7. Conclusions 

In this paper we studied the ES of the spin-| XXZ model in the limit A — > — 1 + 
(ferromagnetic point) for a finite chain with periodic and open boundary conditions. 
To this purpose we used exact diagonalization (with arbitrary precision numerics) and 
Bethe ansatz. We have characterized the ES for all the chain eigenstates within the 
highest spin multiplet (symmetric multiplet) in the limit A — > — 1 + . 
At A = — 1 these states are flat superpositions of all the spin configurations compatible 
with the state magnetization S^. In other terms they show a "mean field" structure 
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with no notion of distance. As a consequence the ES at A = — 1 for each of these 
eigenstates can be calculated exactly and shows few levels. At A > — 1 the "mean field" 
structure is progressively lost as we move away from the ferromagnetic point. 
We summarize our results as follows: 

(i) In the limit A — > — 1 + , most of the ES levels diverge. For each value of the block 
magnetization S\ there is only one non diverging level which reflects components 
of the state wavefunction that in the limit A — > — 1 + are flat superposition of 
all the allowed spin configurations. The divergent behavior is given by £ = 
— alog(A + 1) + (3 with a = 0,2,4,... for periodic boundary conditions and 
a = 0, 4, 8, . . . for open. 

(ii) The multiplicity of an ES level diverging with slope a is given by the sequences 
m per(«) for periodic and m open (a) = m per (a/2) for open boundary conditions. We 
have provided a combinatorial expression [Eq. (12)] for these sequences. 

(iii) We have studied the entanglement eigenfunctions, i.e., the eigenfunctions of the 
reduced density matrix The eigenfunctions of the non diverging ES levels 
mimic the behavior of the state wavefunction, i.e., in the limit A — > — 1 + are flat 
superpositions of all the possible spin configurations. In contrast, the entanglement 
eigenfunctions of the diverging levels do not become flat but display rich structures. 
In particular, in the sector Ma = 1 their leading order in A + 1 can be given 
analytically in terms of the discrete Chebyshev polynomials. 

(iv) Although in the thermodynamic limit in the whole gap less phase at — 1 < A < 1 
the ES distribution is described by (17), for a finite chain significant deviations 
from the CFT prediction and non universal behavior are observed in the vicinity of 
A = —1. Within the range of system sizes considered [La — 8) the ferromagnetic 
point A = — 1 appears to dominate the whole region A < 0. 
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Appendix A. Expansion of exact Bethe ansatz solutions around A = — 1 

In this Appendix we outline the Bethe ansatz expansions used in the Article. 

Appendix A.l. Periodic boundary conditions 

In the region A > — 1 the Bethe equations can be written in the form 

arctan ( ' . \ = —-L + — arctan ( ^ -\ (A.l) 

Vtan (C/2 )/ L L ^—f \(1 - viv m ) t&n( J V ' 

rtvfl 

where M is the number of particles in the chain, (Z = 1,2,..., M) are the so called 
rapidities and ( = arccosA. 1\ are the Bethe numbers and each choice of l\ identifies 
an eigenstate of the XXZ Hamiltonian. In each Sj> sector the lowest energy eigenstate 
(which is the state in the symmetric multiplet with magnetization Sy) is identified by 
the following choice of the Bethe numbers 

/, = _(M + l)/2 + Z / = 1,2,...,M (A.2) 

(Note M = L/2 — Sy.) From the solutions of the Bethe equations the state can be 
constructed as 

(-ire 1 2 «™ m x l7 x 2 ,... t x M ) (A.3) 
v v ' 

a(x 1 ,x 2 ,...,x M ) 

where V runs over all permutations of the set {xi,x 2 , ■ ■ ■ , Xm}, %i being the particle 
positions. We denote with |xi,x 2 , . . . the spin configuration with particles (down 
spins) at sites Xi,x 2 , ■ ■ ■■ The parity of the permutation V is denoted by e-p. We also 
defined ki = 7T + 2 arctan(wz/ tan(^/2)) (mod27r) and the scattering phase k m ) = 
2arctan((t); - u m )/((l - viv m ) tan()). 

The idea of the approach is to solve the Bethe equations (A.l) perturbatively order by 
order in e = A + 1 imposing a solution of the form 

(o) , (i) , 

to obtain a perturbed expression for the state using (A.3). At zeroth order the Bethe 
equations are given by the set of equations 

o i (0) (0) 

„(o) _ V 1 - V l ^ 
' " L ^ V W _ „(o) 

Appendix A. 1.1. First order expansion of the wavefunctions Up to first order in e the 
amplitude (see (A.3)) of the wavefunction with two particles can be written analytically 



a(xi,x 2 ) = Wttt^ — -t 1 + A 2 e - 6 - (x x - x 2 + ^) 2 (A.4) 



L(L-l) 



V 2 
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where 



A, 



L(L-l) 



Xl<X2 



ensures the correct first order normalization. The last term in (A. 4) shows a quadratic 
behavior in the inter particle distance x 2 — #i and linear dependence on e of the 
"curvature" of the wavefunction. Eq. (A. 4) also shows that the asymptotic limit e — > 
the wavefunction becomes an equal weight superposition (no dependence on Xj). These 
features were shown and discussed in Figure 1, where we also compared with exact 
diagonalization data. 

In first order it is possible to go beyond M = 2 and derive the amplitudes for any state 
within the symmetric multiplet: 



a{x 1 ,x 2 , 



M 



l + A M e 



Kj 



Xij + — 



(A.5) 



where X 



M-l 

1/M! n (L - k), X, 

k=0 



X ; 



Xj, and the normalization is 



.4 



11 



(L-l)Z, 



(A.6) 



xi<X2<---<xm K] 

While the M = 2 wavefunction corresponds to an excited state at A > —1, the formula 
above includes as a particular case (M = L/2) the wavefunction of the ground state in 
the gapless phase. One important property of (A.5) is that it is quadratic and symmetric 
in all the interparticle distances Xi—Xj. An example of wavefunction amplitudes obtained 
from (A.5) is shown in Fig. 1 (c). 



Appendix A. 2. Open boundary conditions 

It is possible to treat the open boundary conditions case within the Bethe ansatz 
formalism [24]. The Bethe equations now read 

k ^ Vl) = ITT 1 ' " lTT^' ~ Vl) " IlTi Et $ ^' ~ Vm ^ + $ ^ ( A - 7 ) 

with &(vi,v m ) = 2arctan((t>; — v m )/((l — ViV m )tan()) and k{yi) = ix + 
2 arctan(w; tan(£/2)). The Bethe numbers giving the lowest energy state in each 
magnetization sector (i.e. the states within the symmetric multiplet) are now Ii = 
{1, 2, . . . , M}. The main difference with respect to the periodic boundary conditions is 
the presence of the "reflected" phase shifts $(vi, —v m ). 

This fact has dramatic consequences in the wavefunction and therefore in the ES. The 
general expression of the chain wavefunction has the form 

V 
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where the sum is over all permutations and (arbitrary number of) reflections of the set 
of rapidities {i>i, v 2 , ■ ■ ■ , vm}- The amplitude A(v Vl , v-p 2 , . . . , v-p M ) is not as simple as for 
periodic boundary and we do not show it. As for periodic boundary conditions we can 
expand the Bethe equations in the vicinity of A = — 1. Quite interestingly now we find 
that the lowest order solutions are determined by the equations 



By expanding (A. 8) in A + l we find quite remarkably that the sum over the reflections of 
the rapidities vi — > —V\ cancel at first order in A + 1 the products of the form XiXj. This 
fact leads to the suppression of the first order ES levels for open boundary conditions, 
as we have seen (Figures 2 and 5). We provide as an example of this cancellation the first 
order amplitude a(x\,x 2 ) for the (unnormalized) chain wavefunction in the case with 
two particles, which reads 

a(xx, x 2 ) = 1-4 [L 2 /2 + 3 - L(l/6 + x x + 3x 2 ) + 2x 1 (x 1 + 1) + 2x 2 (x 2 + 1)] (A.10) 

Our numerical results indicate that the same type of cancellation will happen at every 
odd order, but for higher orders this is too cumbersome to perform explicitly. 
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